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Abstract 

We analyze a class of Signal-to-Interference-and-Noise-Ratio (SINR) random 
graphs. These random graphs arise in the modeling packet transmissions in 
wireless networks. In contrast to previous studies on the SINR graphs, we 
consider both a space and a time dimension. The spatial aspect originates from 
the random locations of the network nodes in the Euclidean plane. The time 
aspect stems from the random transmission policy followed by each network 
node and from the time variations of the wireless channel characteristics. The 
combination of these random space and time aspects leads to fluctuations of 
the SINR experienced by the wireless channels, which in turn determine the 
progression of packets in space and time in such a network. 
This paper studies optimal paths in such wireless networks in terms of first 
passage percolation on this random graph. We establish both "positive" and 
"negative" results on the associated time constant. The latter determines the 
asymptotics of the minimum delay required by a packet to progress from a 
source node to a destination node when the Euclidean distance between the 
two tends to infinity. The main negative result states that this time constant 
is infinite on the random graph associated with a Poisson point process under 
natural assumptions on the wireless channels. The main positive result states 
that when adding a periodic node infrastructure of arbitrarily small intensity 
to the Poisson point process, the time constant is positive and finite. 
Keywords: Poisson point process, random graph, first passage percolation, shot- 
noise process, SINR 
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1. Introduction 

There is a rich literature on random graphs generated over a random point process. 
These graphs are often motivated by physical, biological or social networks. Many 
interesting large scale properties of these networks related to connectivity have been 
studied in terms of the percolation of the associated graphs. An early example of such 
a study can be found in [13] where the connectivity of large networks was defined as 
the supercritical phase in what is today called the continuum (Boolean) percolation 
model. More recently, a random SINR graph model for wireless networks was studied 
with the same perspective in [10, 11]. 

The routing, and more precisely, the speed of delivery of information in networks 
is another example of problems, which motivated the study of the random graphs. 
The main object in this context is the evaluation of the so called time constant, which 
gives the asymptotic behavior of the number of edges (hops) in the paths (optimal or 
produced by some particular routing protocol) joining two given nodes in function of 
the (Euclidean) distance between these nodes, when this distance tends to infinity. In 
the case of a shortest (in terms of the number of hops) path, this problem is usually 
called the first passage percolation problem and was originally stated by Broadbent and 
Hammersley in [8] to study the spread of fluid in a porous medium. More recently, 
in [20, 7], such time constants were studied on so called small word graphs, motivated 
by routing in certain social networks, where any two given nodes are joined by an 
edge independently with a probability that decays as some power function with the 
Euclidean distance between them. The complete graph on a Poisson p.p. with "nearest 
neighbor" routing policy was studied in this context in [6] . The first passage percolation 
problem on the Poisson-Dclaunay graph was considered in [24, 21]. In the case of graphs 
whose edges are marked by some weights, one can extend the notion of time constant 
by studying the sum of the edge weights. First passage percolation on the complete 
Poisson p.p. graph, with weights proportional to some power of the distance between 
the nodes was studied in [14]. 

The present paper focuses on the speed of delivery of information in SINR graphs. In 
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contrast to previous studies of this subject, in particular to [10, 11], we consider graphs 
with space-time vertexes. This new model is motivated by multihop routing protocols 
used in wireless ad-hoc networks. In this framework, the random point process on 
the plane describes the locations of the users of an ad-hoc network and the discrete 
time dimension corresponds to successive time slots in which these nodes exchange 
information (here packets). As in [2] we assume the spatial Aloha policy to decide 
which node transmit at a given time slot. We also assume some space-time fading 
model (already used e.g. in [5]) to describe the variability of the wireless channel 
conditions (see e.g. [ -'■■]). In this space-time SINR graph, a directed edge represents 
the feasibility of the wireless transmission between two given network nodes at a given 
time. More precisely, the direct transmission of a packet is succeeds between two 
nodes in a given time slot if the ratio of the power of the signal between these nodes 
to the interference and noise at the receiver is larger than a threshold at this time 
slot. This definition has an information theoretic basis (see e.g. [23]). It is rigorously 
defined below using some power path-loss model and an associated shot-noise model 
representing the interference. 

We study various problems on this random graph including the law of its in- and 
out-degree, the number of paths originating from (or terminating at) a typical node 
or its connectedness. The most important results bear on the first passage percolation 
problem in this graph. In the case of Poisson p.p. for the node locations, we show 
that the time constant is infinite. We then show that when adding a periodic node 
infrastructure of arbitrarily small intensity to the Poisson point process, the time 
constant is positive and finite. These results lead to bounds on the delays in ad-hoc 
networks which hold for all routing algorithms. This subject, or more generally, the 
question of the speed of the delivery of information in large wireless ad-hoc networks 
currently receives a lot of attention in the engineering literature see e.g. [12, 15]. 

The paper is organized as follows. In Section 2 we introduce the space-time SINR 
graph model. The results are presented in Section 3, Most of the proofs are deferred to 
Section 4. Some implications on routing in ad-hoc networks are presented in Section 5. 
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2. The Model 

2.1. Probabilistic Assumptions 

Throughout the paper we consider a simple, stationary, independently marked (i.m.) 
point process (p.p.) $ = {(Xi, e,;, F;, W^)} with finite, positive intensity A on K 2 . In 
this model, 

• $ = {Xi} denotes the locations of the network nodes on the plane M 2 . The 
following three cases regarding the distribution of $ will be considered: 

General p.p.: $ is a general (stationary, non-null, with finite intensity) p.p., 
Poisson p.p.: <f> is a Poisson p.p., 

Poisson+Grid p.p.: <& = &m + <&g is the superposition of two independent 
p.p.s; where &m denotes a stationary Poisson p.p. with finite, non-null 
intensity Xm and <&g = sZ 2 + Ug a stationary, periodic p.p., whose nodes 
constitute a square grid with edge length s, randomly shifted by the vector 
Ug that is uniformly distributed in [0, s] 2 (this makes $g stationary). Note 
that the intensity of $g is = 1/s 2 - 

• e, = {e,(n)} 6ne z, where % = {. . . , —I, 0, 1, . . .} denotes integers; the variables 
{ei(n) : are i.i.d. (in n and i) Bernoulli random variables (r.v.s) with 
P{ e = 1 } = 1 — P{ e = } = p, where e denotes the generic r.v. for this family. 
We always assume < p < 1. The variable ei(n) represents the medium access 
indicator of node Xi at time n; it says whether the node transmits (case e = 1) 
or not at time n. 

• W t = {Wi(7i)}„; {Wi(n) : i,n} is a family of non-negative i.i.d. r.v.s with some 
arbitrary distribution. The variable Wi(n) represents the power of the thermal 
noise at node X t and at time n. Let W denote the generic r.v. for this family. 

• Fj = {Fi,j(n)}j t n\ {Fij(n) : i,j,n} is a family of non-negative i.i.d. r.v.s. The 
variable Fij(n) represents the quality of the radio channel (also called fading) 
from node Xi £ $ to node Xj G $ at time n. The following two cases regarding 
the distribution of F (denoting the generic random variable for this family) will 
be considered: 
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General fading: when F has some arbitrary distribution with finite mean. 
Exponential fading: when F has exponential distribution with mean 1/yU. (*) 

To complete the probabilistic description of the model we assume that, given <£, the 
random elements {e^};, {W^}, and {F^}.; are independent. For more on this framework, 
which is classical, see e.g. or [3, 4, 5]. 

Our stationary i.m. p.p. $ is considered on some probability space with probability 
P. We will denote by P° the Palm probability with respect to see [9, Ch.13]. Recall 
that it can be interpreted as the conditional probability given $ has a point at the 
origin of the plane. We will denote this point (considered under P°) by Xq and call 
it the typical node. Under P° $ is also an i.m. p.p. with marks distributed as in the 
original law. Moreover, in the case of Poisson p.p.s, the distribution of <& under P° is 
equal to the distribution of $ U {X = 0} under the stationary probability P (cf. the 
Slivnyak-Mecke Theorem [9, p.281]). 

2.2. SINR Marks 

Given the i.m. p. p. $ described above, we construct another family of random 
variables {SINRjj(n) : n}, which will be interpreted as the SINR observed in the 
channel form Xi £ $ to Xj £ $ at time n. These variables, which have an information 
theoretic background, will be used to assess the success of transmissions. For defining 
these variables, we give ourselves some non-decreasing function I : K+ = {t : t > 0} — ► 
R + that we call the path-loss function. A special example considered in this paper 
(and commonly accepted in the wireless communication context) is 

l(r) = (Arf with some A > and (3 > 2 . (1) 

Denote by $ 1 (n) = {Xi : e,(n) = 1} the point process of transmitters at time slot n 
and by <&°(n) = {Xi : ei(n) = 0} that of (potential) receivers. For a given pair receiver 
Xj £ $° and transmitter Xi £ $ 1 (n), we will assume that Xj receives a signal from Xi 
with power Fij(n)/l(\Xj — Xi\) at time n. Node Xj also receives signals from other 

(*) In wireless signal propagation models, the exponential distribution appears naturally as the 
square power of the norm of a complex random variable, whose real and imaginary components are 
i.i.d. Gaussian. I this case one often speaks about Rayleigh fading model because the norm (absolute 
value) of such a complex random variable is Rayleigh-distributed; see e.g. [23, p. 50 and 501]. 
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transmitters X k € 3? (n), Xf. ^ Xi at time n. The total received power is equal to 
Let also 

Ij(n)= E ^HA(|X fe -X,|). 

X fe g*i(")\{^} 

Both Iij(n) and ij(n) are shot-noise r.v.s generated by $ 1 (n), the fading marks and 
the path-loss function. The are infinite sums of non-negative r.v.s. In order to check 
whether these r.v.s are a.s. finite, one can use the Campbell-Little-Mecke formula 
(Campbell for short; cf. [9, Prop. 13.3.II] ) , which implies that 



E F kfi (n)/l(\X k \)]= P B[F] [ l/l(\x\)M [2] (dx) , (2) 



where M2O) is the reduced second order moment measure of $ (cf [9, p. 238]). In 
what follows, we will always tacitly assume that & are such that the integral in 
the right-hand side of (2) is finite for some e > 0, which implies that Io(n) is almost 
surely (a.s.) finite under P° for all n as well as all Ij(n), h,j{n) under P. If $ is the 
homogeneous Poisson p.p., we have M[ 2 ](dx) = Ada; and it is easy see that the we have 
finiteness for /(•) given by (1) for all e > 0. It is also relatively easy to see that it holds 
for the Poisson+Grid p.p. $ = $m + $g- 

The SINR at the receiver Xj G &°(n) with respect to transmitter Xi G ^(n), at 
time n is defined as 

» jH =^ta. <3) 

Wj{n)+Ii,j[n) 

2.3. Space-Time SINR Graph 

Let 

l(SINRi,j > T) if e l {n) = 1, ej (n) = 0, i ^ j, 

kj( n ) = \l if*=j, ( 4 ) 

otherwise, 



Optimal Paths on the Space-Time SINK Random Graphs 



7 



where T > is a threshold assumed to be some given constant throughout the paper. 
We define the space-time SINR graph G as the directed graph with the set of vertexes 
$ x Z and a directed edge from (X,-, n) to (Xj, n + 1) if Sij (n) = 1. 

Let us stress an important convention in our terminology. By network node, or 
point, we understand a point of <f>. A (graph) vertex is an clement of $ x Z; i.e. it 
represents some network node at some time. The existence of a graph edge is to be 
interpreted as the possibility of a successful communication between two network nodes 
(those involved in the edge) at time n. This can be rephrased as follows. Suppose that 
at time n the network node Xi has a packet (containing some information) . Then the 
set of graph neighbors of the vertex (Xj , n) describes all the nodes that can decode this 
packet at time n + 1 . Thus any path on the graph G represents some possible route of 
the packet in space and time. 

3. Results 

In this section we present our results on G. 

3.1. Existence of Paths 

All the results of this section arc obtained under the general p.p. and fading 
assumptions of Section 2, under the assumption that the finiteness of the expression 
in (2) is granted. 

Note first that G has no isolated nodes in the usual sense. Indeed, we have always 
(Xi, n) connected to (X^, n + 1). We will consider directed paths on G and call them 
paths for short. Note that these paths are self-avoiding due to the fact that there are 
no loops in the time dimension. 

Denote by Ti.° ut ' k (n) the number of paths of length k (i.e. with k edges) originating 
from (A,j,n). Similarly, denote by 7^ n ' (n) the number of such path terminating at 
(Xi,n). In particular H° ut {n) = H" ut,1 (n) and H\ n (n) = H 1 ™' (n) are respectively, the 
out- and in-degree of the node (Xi,n). 

Lemma 1. For a general p.p. $ and a general fading model, the in-degree H\ n of any 
node of G is bounded from above by the constant £ = 1/T + 2. 

Proof. Assume there is an edge to node (Xj , n) from nodes [X^ , n— 1), . . . , (Xi k , n — 
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1), for some k > 1 and i p ^ j (p = 1, . . . , /c). Then for all such p 

z(|x ip - - TTr \Y i(|*, 4 - ; ' 

When summing up all these inequalities, one gets that Tk < 1 + T, that is k < 1/T+l. 
Considering the edge from (X,, n — 1) to (X,;, n), the in-degree of any node is bounded 
from above by £ = l/T + 2. □ 

Let 

h out,k = E 0[ W ^,fc( n )] = E°[W™*' fe (0)] 

and 

fi in -' = E°[H™' fe (n)] = E°[74"' fc (0)] 

be the expected numbers of paths of length k originating or terminating at the typical 
node, respectively. In particular h out — /j™*- 1 and h m = h" hl are the mean out- and 
in-degrec of the typical node, respectively 

Lemma 2. For a general p.p. $ and a general fading model 

h in,k = h out,k ^ 

Proof. We use the mass transport principle to get that E° [Ti.Q Ut ' k (0)] = E° [Hg ra ' fc (0)] , 
which implies the desired result. Indeed, Campbell's formula and stationarity give 



Xh ouLk = X ( E°[H™*' fc (0)]da; 

J[0,1) 2 



e[ K ut - k (o)} 

Xie<s>n[Qj) 2 



Y E [ Y Y # of paths from ( x * ' °) to ' k ^> 

vez Xie[o,i) 2 Xje[o,i) 2 +u 

Y E [ Y Y # of paths from ( Xi , 0) to (JiTj, fc) 



= A 



i[o,i) 2 



where # denotes the cardinality. This completes the proof. □ 
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Here are immediate consequences of the two above lemmas. 
Corollary 3.1. Under the assumptions of Lemma 1: 

• G is locally finite (both on in- and out-degrees of all nodes are P-a.s. finite). 

• n^ k (n) < £ fc P-a.s for all i, n, k. 

• h m > k = h out ' k < £ fe for all k. 

For all Xi, X, G $ and n G Z, we denote by we will call local delay from Xi to Xj 
at time n the quantity 

Li.j{n) = infjfc > n : 8ij(k) = 1} 

with the usual convention that inf = oo. Note that Lij(n) is the length (number 
of edges) of the shortest path (with the smallest number of edges) from (X, , n) to 
{Xj} x Z among the paths contained in the subgraph G n {X i: Xj} x Z of G, which is 
of the form 

((X u n), (Xi,n + 1)), . . . , {{Xi,n + L id (n) - 1), (X u n + L itj {n))), 
((Xi,n + Li^n)), (X j>n + L itj {n) + 1)). 

Our next result gives a condition for the local delays to be a.s. finite. 

Lemma 3. Assume a general p.p. $ and a general fading model with F having 
unbounded support (P{ F > s} > for all < s < oo). Then, given all local 
delays Li.j{n) are P-a.s. finite geometric random variables. 

Proof. Due to our assumption on the independence of marks in successive time slots, 
given <3>, the variables {8ij(n) : n G Z} are (i.i.d.) Bernoulli r.v. and thus Lij(n) is 
geometric r.v. It remains to show that P{ Sij(0) = 1 1 <J> } := 7r,; j($) > for P-almost 
all $. For this, note that 

TT^*) =p(l - p )p{f ij (o) > «(|X i -^|)(w j (0)+/ iJ (0)) } . 

Under our general assumptions (including finiteness of the expression in (2)) Ii.j(O) is 
a finite random variable P-a.s. The result follows from the assumption that < p < 1 
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and the fact that Fij(0) is independent of lij(0), Wij(0) and has infinite support. □ 

The next result directly follows from Lemma 3. 

Corollary 3.2. Under the assumptions of Lemma 3, G is P-a.s. connected in the 
following weak sense: for all Xi, Xj G <£> and all neZ, there exists a path from (Xi, n) 
to the set {(Xj,n + l) : I G N}, where N = {1,2,.. .}. 

We denote by Li(n) = infjyj L^j the length of a shortest directed path from (Xi, n) 
to ({$ \ Xi}) x Z. We will call Lj(n) the exit delay from Xi at time n. Finally, we 
denote by Pij(n) the length of a shortest path of G from (Xj, n) to {Xj} x Z. We call 
Pij(n) the deZa?/ from Xi to X,- at time n. Obviously for i =/= j we have 

Li{n) <Pi,j(n) <L id (n) (6) 

and thus it follows immediately from Lemma 3 that all the three collections of delays 
finite r.v.s P-a.s. 

3.2. Optimal Paths — Poisson p.p. Case 

We have seen in the previous section that under very general assumptions, all the 
delays arc P-a.s. finite random variables. In this section we show that under some 
natural assumptions (such as Poisson p.p. and exponential fading), the averaging 
over $ may lead to infinite mean values. This averaging is expressed in terms of the 
expectation for the typical node under the Palm probability. The proofs of the results 
stated in what follows are given in Section 4.1. 

Denote t = E°[L (n)] = E [L o (0)]. 

Proposition 3.1. Assume $ to be a Poisson p.p., F to be exponential and the noise 
W to be bounded away from 0: P{ W > w} = 1 for some w > 0. Let the path-loss 
function be given by (1). Then P°{io(0) > q} > l/q for q large enough. 

Corollary 3.3. Under the assumptions of Proposition 3.1, we have: 

• The mean exit delay from the typical node is infinite; I = oo. 

• In any given subset of plane with positive Lebesgue measure, at a given time, 
the expected number of points o/$ which have exit delays larger than q decreases 
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not faster than 1/q asymptotically for large q. 

The fact that the mean exit delay from the typical point is infinite (t = oo) seems to 
be a consequence of the potential existence of arbitrarily large "voids" (disks without 
points of <f>) around this point. Indeed, when conditioning on the existence of another 
point in the configuration $, one obtains finite mean local delays. This will be shown 
in Proposition 3.2 below. 

Before stating it we need to formalise the notion of existence of two given points 
X, Y £ M 2 of For this, we consider $ under the two-fold Palm probability P X ' Y . 
Since our results on the matter bear only on the Poisson p.p. case, we can assume (by 
Slivnyak's Theorem) the following version of the Palm probability of the Poisson p.p. 
i>: 



Moreover, under P ' , the marked Poisson p.p. <!> is obtained by an independent 
marking of the points of $ U {X, Y} according to the original distribution of marks. 
Slightly abusing the notation, we denote by Lx,Y{n) the local delay from X to Y 
at time n when considered under P X,Y . Similar convention will be adopted in the 
notation of other types of delays under the Palm probabilities or P X,Y . 

Proposition 3.2. Assume $ to be a Poisson p.p., F to be exponential and the noise 
W to have a general distribution. Then for all X,Y £ R 2 , the mean local delay from 
X to Y is finite given the existence of these two points in (f>. More precisely, 




(7) 



E x < Y [L x , Y (0)] <oo. 



(8) 



The next result follows immediately from (6). 



Corollary 3.4. Under the assumptions of Proposition 3.2, 



E x < y [L x (0)] < E x > Y [P x ,y(0)] < E X > Y [L XX (0)] < oo. 



(9) 



The following result is our main "negative" result concerning G in the Poisson p.p. 
case: 
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Proposition 3.3. Under the assumptions of Proposition 3.1, we have 

| Ace \X-Y\ ^ ■ (10) 



In other words, the expected shortest delay necessary to send a packet between two 
given points of the Poisson p.p. grows faster than the Euclidean distance between these 
two points. 

3.3. Filling in Poisson Voids 

In this section we show that adding an independent periodic pattern of points to the 
Poisson p.p. allows one to get a linear scaling of the shortest path delay with Euclidean 
distance. In order to prove the existence and finiteness of the associated time constant, 
we adopt a slightly different approach to the notion of paths on G, which will allow 
us to exploit a subadditive ergodic theorem. The proofs of the results stated in what 
follows are given in Section 4.2. 

For ieE 2 , let X(x) be the point of <E> which is closest to x. The point X{x) G $ 
is a.s. well defined for all given x £ R 2 since $ is assumed simple and stationary p.p. 
For all x, y € M 2 , define P(x, y, n) = Px(x),x(y)( n ) to be the length of a shortest path 
of G from vertex (X(s), n) to the set {(X(y), n + 1), I G N}. We will call P(x, y, n) the 
delay from x to y at time n. For all triples of points x, y, z G R 2 , we have 

P(x, z, n) < P(x, y, n) + p(y, z,n + P(x, y, n)^j . (11) 

Let 

p(x,l/ I *)=E[P(x J »,0)|#]. (12) 

Using the strong Markov property, we get that, conditionally on $, the law of P(y, z, n+ 
P(x, y, n)) is the same as that of P(y, z, n). Then, the last relation and (11) give 

p(x,z,$)<p(x,y,$)+p(y,z,$). (13) 



We are now in a position to use the subadditive ergodic theorem to show the 
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existence of the time constant 

K d = lim , 

t—>oo t 

where may depend on the unit vector d € I 2 representing the direction in which 
the delay is measured. Here is the main result of this section. 

Proposition 3.4. Consider the Poisson+Grid p.p. defined in Section 2.1 with expo- 
nential jading F and with the path-loss function be given by (1). Then, for all unit 
vectors d £ M. 2 , the non-negative limit Kd exists and is P-a.s. finite. The convergence 
also holds in L\ . 

Notice that k& is not a constant. Indeed, the superposition of the p.p.s $ = $m 
and is ergodic but not mixing due to the fact that $ is a (stationary) grid. For 
d parallel to say the horizontal axis of the grid ^g, the limit Kd will depend on the 
distance from the line {id : t € M} to the nearest parallel (horizontal) line of the grid 
$G, i.e. on the shift Uq of the grid. Here is a more precise formulation of the result. 

Proposition 3.5. Under the assumptions of Proposition 3.4, the limit Kd = Kd(Ua) is 
measurable w.r.t. the shift JJq of the grid p.p. $g o,nd does not depend on the Poisson 
component $m of the p.p. <f>. Moreover, the set of vectors d in the unit sphere for 
which KdiUc) is not P-a.s. a constant is at most countable. 

The last result on this case is: 

Proposition 3.6. Under the assumptions of Proposition 3.4, suppose that W is con- 
stant and strictly positive. Then E[«d] > 0. 

Finally let us remark that the method used in this section cannot be used in the 
case of the Poisson p.p. (without the addition of the grid point process). The main 
problem is the lack of integrability of p(x, y, $) as stated in the following result. Note 
however, that this does not imply immediately that Kd = oo. 

Corollary 3.5. Under the assumptions of Proposition 3.1 E[p(a;, y, $)] = oo for all x 
and y in Mr. 
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4. Proofs 

Consider the shortest path from (Xi,ri) to ($\ {^i}) x Z. Let %(n) be the number 
of edges (Xi, k), (X,-, k + 1) in this path such that ej(fe) = 1. These variables are the 
number of trials before the first exit form Xi at time n. Obviously 



%{n) < Li(n) . (14) 

We will also consider an auxiliary graph G, called the (space-time) Signal to Noise 
Ratio (SNR) graph, defined exactly in the same manner as the SINR graph G except 
that the variables SINRij(n) defined in (15) are replaced by the variables 

SNRjj(n) = fkWZfcM. (15) 

Note that this modification consists in suppressing the interference term I^ j (n) in the 
SINR condition in (4). The edges of G form a subset of the edges of G (both graph 
share the same vertexes) , which will be denoted by 

GcG. (16) 



In what follows we will denote the delays, local delays, exit delays and numbers of trials 
related to G by Pij(n), Lij(n), Li(n) and %(n), respectively. The inclusion G C G 
implies immediately that Pij(n) < Pij(n) and the same inequalities hold for the three 
other families of variables mentioned above. 

4.1. Proofs of Results of Section 3.2 

Proof, (of Proposition 3.1) The inclusion (16) and the inequality (14) yield 

%(n) < %{n) < Li(n) , 

which holds for all i, n. The results follow from the above inequalities and the next 
lemma. □ 

Lemma 4. Under the assumptions of Proposition 3.1, P {7q(0) > q } > 1/q for q 
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large enough. 



Proof. Under P°, denote by Tt the kth time slot in {0, 1, . . .}, such that eo(k) = 1. 
For all q > we have 

P°{f o (0) > = P o {Vo< fc < g V o# x ie *^(T fc ) = 0|$} 

= P°{ V < fc < g V 0#M e^n-) = 1 or SNR ^(r fc ) < T | $ } 

and by the conditional independence of marks given $ 

P°{f o (0)>g|$} = [] (P+(1-P)P{^<^(|X 4 |)^}) 9 

= exp{, J2 log(p+(l-p)(l- e -™ l)H/ ))}, 

where i 7 ", W are independent generic random variables representing fading and thermal 
noise, independent of 3>, i* 1 is exponential with mean l//t, Using the Laplace functional 
formula for 4> and the assumption that W > w a.s. we have 

P°{%(Q)>q} > cxp^-27rA^ (l - (l - (1 - p)e- w ^ T Y^ v dvj 



= cxp t^-nX y (1 - (1 - /(„))«) d«J , (17) 

where 

/(u) := (1 - p) exp(-Kv p/2 ) and A' = iu/xTA^ . 

In what follows we will show that the expression in (17) is not smaller than 1/q for q 
large enough. To this regard denote by v q the unique solution of f(v) = |. We have 



It is clear that f(v) tends to when v tends to infinity and that v q tends to infinity 
as q tends to infinity. Therefore, there exists a constant Q = Q(fi,w, A,T) < oo such 
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that for all q > Q and for all v > v g , 

(1 - /(«)) > exp(-/(«)). 

Hence, for all q > Q : 

pOO 

(1 - (1 - f(v)) q ) dv < v q + (1 - (1 - f(v)) q ) dv 



v>0 

< v q + I (1 - exp(-qf{v)) dv 



v=v q 
oo 



< V q + / qf(v) dv 



/•OO 

= V g + qf(u + v q ) du. 

Ju=0 

The third inequality follows from the fact that 1 — exp(— x) < x. Using now the fact 
that (m + Vq) 13 / 2 > u + Vq^ 2 (for q large enough, say again q > Q) we get that 



qf{u + v q )du = / q(l - p) cxp(-K(u + v q f /2 j) du 

Ju=0 

f°° 1 

< / q(l - p) cxp( — I\u - Kv^ 2 ) du = — 

Ju=Q K 



u=0 Ju=0 

oo 



Q /2 

since (1 — p) exp(—Kv q ) = 1/q. Hence for q > Q 



(l-(l -/(«))*) dv<V q + ^. 
v>0 " 



Also it is not difficult to see that (3 > 2 implies 



V ^^X-K (18) 



(19) 



for q large enough. This implies for q large enough, say again q > Q, 

cxp (-nxj (1 - (1 - f(v)) q ) dv^j > exp(-7rX(v q + 1/A')) > ^ , 

which completes the proof. □ 
Proof, (of Proposition 3.2). Assume without loss of generality Y = and |A| = r. 
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Under P, consider the p.p. $ U {X, 0} and its independent marking. Given the r.v. 
Lx,o{0) associated with the independently marked p.p. <& U {^,0} has a geometric 
distribution with parameter 

7r X) o($) = p(l - p) Pr{ F > l(r)(W + 1))}, 

where F, W, I arc independent r.v.s. F, W are generic fading and noise variables and 
I = Ylx e<s> e i(0)-^,o(0)/£(|^Q|)- Using the exponential distribution of F and the 
independence, we obtain 

n Xfi (<P) = E[e~^ TW ] Efe-'W | <f>] . 

The mean of the geometric r.v. is known to be E X '°[Lx,o(0) | $] = l/irx,o(&)- By 
unconditioning with respect to <f>. one obtains 



E^ u [L x , o (0)] = 



1 



-E 



1 



£ w (pl(r)T) LE[e-"'M 7V |$] 



The first factor in the above expression is obviously finite. In what follows we will 
evaluate the second one. 

By the conditional independence of marks and denoting by C £ f{-) is the Laplace 
transform of eF, where e, F are independent generic variables for ei(0) and -Fi.o(O) we 
have 



(E[e-^' T/ |$]]) = [e 



exp(-/iZ(r)T £ e*(0)F <lO (n)/J(|^i|) 



exp( J] log£ e F (juT/(r)/Z(|Xi|)) J 



Note that C b f(£) = 1 — p + p£f{0 = 1 — P + P/V (/■* + £)• Using this and the Laplace 
functional formula for $, (cf. [9, Eq. 9.4.17]) we obtain 



E 



1 



E [ e -Mi(r)T/ | $] 



= exp< 2 7175 A 



i/H(r) 



o J(«) + (l-p)Ti(r) 



dc 



(cf. (2) Using now the fact that for the Poisson p.p., M\2](dx) — Ada;), it is now easy 
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to see that for any path-loss function satisfying v/l(v) dv < oo, the integral in the 
exponent of the last expression is finite. This completes the proof. □ 

Proof, (of Proposition 3.3). Using the the inclusion (16), inequality (14) and the 
left-hand side of (6) and we have 

fi(n) <%{n)< Li(n) < P w (n) . 

Thus, it is enough to show 

hm — — — — = oo . 

|x-y|->-oo 

Without loss of generality assume X = and \Y\ = r. Using the same arguments as in 
the proof of Lemma 4 and the representation (7) of the Palm probability with respect 
to Poisson p.p., we obtain 

P°< y {f (0)> g |<I>} 

> II (p+(1-p)P{F <Tl(\X l \)W}) q (p+(l-p)P{F<Tl(\Y\)W} 



> cxp(-n\J^ (1 - (1 - f(v)) q ) dv) a{r) q , 



where a(r) = 1 — (1 — p)e wfJ/A TrP . Using (19), which holds for large q, more precisely 
q > Q = Q(n, u>, A, T), we obtain 

~ r h q 

It is now easy to see that 

1 v-^ a ( r ) q 
hm - } = oo. 

q>Q H 

□ 



4.2. Proofs of Results of Section 3.3 

Denote by B X (R) the ball centered at x £ R 2 of radius R. Similarly as for the delays, 
we extend the definition of the local delays to arbitrary pairs of points x, y £ U 2 by 
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taking L(x,y,n) = Lx( x ),x(y){ n )- We first establish the following technical result: 

Lemma 5. Under the assumptions of Proposition 3.4 let Xi,Xj G ^CiBo(R) for some 
R > 0, where 4> = <&m + ^g 3 - Then the conditional expectation of the local delay Ljj(O) 
given $ satisfies 

E[ii,i(0)|$] 

'T\Xi-XA^- 



{ E ^40} 



pil-p)C w {T^\Xi - ) GXP 1 $ 3^ )fe ^ 
1 

< 



p{l-p)C w {Tp(A2R)P) 

x e -491og(l-p) + (2i?.)' 5 P TC(s,f3) 



Xe -#M(B (2R))log(l-p) 



x exp 

^ fe e*A/,l^l>2fl 



(|X fc |-i?)^+T(2i^ 



where C(s,f3) < oo is some constant (which depends on s and [3 but not on <f>), F' is 
an exponential random variable of mean 1 and C e F'{-) is the Laplace transform of eF' . 

Proof. We first prove the equality in (20). When using the independence assump- 
tions, we have 

P{L itj (0) >m\ $} 

= P{V^ =1 (e i (n) = lor 

e,{n) = and e^F^n) < Tl(\X t - X,|)(V^(n) + h.^n))) } 

= J I P + (1 - p) ( 1 - P + p(l - £w{T»AP\x - yf) 



n c -f' 

<5>3X k ,k^i,j 



T\X t -Xjf 



The result then follows from the evaluation of 



E[L^(0) | <E>] = ^P{L w (0)>m | $] . 

m=0 

The bound |Xj — X,| < 2i? used in the Laplace transform of W leads to the first 
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factor of the upper bound. We now factorize the exponential function in (20) as the 
product of three exponential functions 



a : = expj- ^ j, 

(3 : = expj- }> 

$M3X k ,k^i,j\X k \<2R 

7 := exp{- }■ 



$M3X k ,\X k \>2R 

Next we prove that the last three exponentials are upper-bounded by (a), (b) and (c) 
in (20), respectively. 



(a) We use \Xi — Xj\ < 2R and Jensen's inequality to get 

log £ eF , — -^s- > log CeF' ' 



x 3 -x k \p ) - \\Xi-X k \» 

T(2R)PE[eF'} 



> 



\Xj-X k \P 
-pT(2Rf\X 3 - X k \~ p 



We now prove that 



]T \Xi-Xu\-P <C{s,p), 

®G e 3X k :\Xj-Xk\>3V2s 

for some constant C(s, (3). This follows from an upper-bounding of the value of 
\Xj - X k \~ by the value of the integral 1/s 2 J(\Xj - x\ - \/2s)~' 9 dx over the 
square with corner points X k , X k + (s,0), X k + (0, s) and X k + (s, s). In this 
way we obtain 



E M-X^ < ~ (\Xj-x\-V2s)-? dx 

s J\x-X j \>2^/2s 

2tt f°° t + V2s 



. p dt=:C(s,p) <oo. 
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Combining this and what precedes, we get that 

'T\Xj-Xtf* 



cxp j E log£ eJ v( ) Ucxp(T(2i?)' 3 C( S ,/3)). 



X fc 6$ Gs ,l^j-^ fc |>2v^ 



We also have 



logCeF' L _ x , J >log£ e F'(°°) =log(l-p), 

for all Xfc g $g s and in particular for \Xj — Xk\ < 3\/2s. Hence we obtain 



cxp{- E (•")}< 



e -491og(l-p)+T(2fl)' J C( S ,/3) 



X k et>c 



where 49 upper-bounds the number of points Xk G $g s such that |Xj — Xk\ < 
3V2s. 

(b) Using the bound \Xj— Xi\ < 2R and the inequality log C e F' (0 > log£ e F'(oo) = 
log(l — p), we obtain 

CX P {- E (...)} < e -**(*o(2K)) log(l-p) _ 

(c) Using the bounds \Xj-Xi\ < 2R and \Xj-X k \ > \X k \-R (the latter follows 
from the triangle inequality) and the expression C e F'(0 = 1 — p + jfe, we obtain 



cxp{- £ (...)} 

<S> M 3X k ,\X k \>2R 



This completes the proof. 



□ 



We can now prove the following auxiliary result. 
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Lemma 6. Under the assumptions of Proposition 3.4 for all points x,y of] 



E 



sup p(xi,y\,$) 



< 00 , 



where the supremum is taken over Xi,y\ belonging to the interval [x,y] C M 2 . 

Proof. Without loss of generality, we assume that (x + y)/2 = is the origin of 
the plane. Let B = Bq(R) be the ball centered at and of radius R such that no 
modification of the points in the complement of B modifies X(z) for any z G [x,y] 
(recall that X(z) is the point of $ which is the closest from z). Since <& = $m + $g s > 
with the square lattice p.p. with intensity 1/s 2 , it suffices to take R = \u — v\/2 + 
\/2s. Let B' = Bq(2R). By the above choice of B and the inequality (6) we have for 
all xt,yi G [x,y] 

P(xi,yi,0)< J2 

Xi,x 3 e<i-nB 

and consequently 



sup p{x u y u $)< J2 E [ L «(°)I $ ]- 
Using the result of Lemma 20 we obtain 



sup |p* (£1,1/1,$) | 

xi,yiE[x,y] 

e -49 log(l-p) + (2fl)' 3 pTC(s,/3) 

" p{\-p)C w {T^A{2RY) 

v / \^ ^ J^ p{\X k \-Rf ^ 

x($ M (B) + 7 r(i? + \^ s ) 2 / s 2 )e-*«( B ')iog(i- P ) ; 

where n(R+V2s) 2 /s 2 is an upper bound of the number of points of $g s in B. The first 
factor in the above upper bound is deterministic. The two other factors are random 
and independent due to the independence property of the Poisson p.p. The finitcness 
of the expectation of the last expression follows from the finitcness of the exponential 
moments (of any order) of the Poisson random variable <f>j\/ (£?'). For the expectation 
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of the second (exponential) factor, we use the known form of the Laplace transform of 
the Poisson SN to obtain the following expression 



E 



exp{- £(...)}] =cxp|2^A M vP + [l _ p)Tm ^ + R) eh} 



< CO . 



□ 



Proof, {of Proposition 3.4 The existence and finiteness of the limit Kd follows from 
the subadditivity (13) and Lemma 6 by the continuous-parameter sub-additive ergodic 
theorem (see [19, Theorem 4]). □ 

Proof, (of Proposition 3.5) First, we prove the second statement; i.e., that Kd is 
constant for all d in the unit sphere off some countable subset. Note that the point 
process $ is ergodic as the independent superposition of mixing Poisson p.p. $m 
and ergodic grid process <&g- This can be easily proved using e.g. the respective 
characterisations of above properties by means of Laplace transforms of p.p. (see [9, 
Prop. 12. 3. VI]). From the ergodicity of <& we cannot conclude the desired property for 
any vector d since the limit Kd — Kd($) is n °t necessarily invariant with respect to 
translations of $ by any vector x € IR 2 but only x = ad for any scalar a g R. The 
announced result follows from [22, Th. 1]. 

For the first statement, consider a product space with on which two independent 
p.p.s ($a/, <J>g) are defined. Fix some vector d and define the operator T = T\ x 
on this product space as the product of two operators, which correspond to the shift 
in the direction d, of $m and 4>g respectively. The cr-field invariant with respect to 
T is the product of the respective cr-fields invariant with respect to T\ and T%. The 
latter is trivial since <&m is mixing (as a Poisson p.p.). Consequently every function 
of $g) that is invariant with respect to the shift in the direction d of its first 
argument (3>m) is a.s. constant. This concludes the proof that Kd is constant in <3?m 
and thus depends only on Ug- 

Proof, (of Proposition 3.6) For a given path a = {(Xq,uq), (Xi,no+l), ■ ■ ■ , (Xk,no+ 
k)} on G denote by \a\ = S»=i l-^i ~ the Euclidean length of the projection of 

a on R 2 ; let us call it Euclidean length of a for short and recall that the (graph) length 
of a is equal to fc. For fixed e > and all n > 1 denote by II(n) = n e (n) the event that 
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there exists a path on G starting at (X(0),0) that has (graph) length n and Euclidean 
length larger than n/e. 

Assume E[«a] = 0- We show first that this implies that for any e > 0, P°-a.s. the 
event II e (n) holds for infinitely many n 

p°{n |Jn e (fc)}=i. (21) 

n> 1 k>n 

Indeed. E[/«d] = implies = P-a.s. and by Palm-Matthcs definition of the Palm 
probability P°-a.s. as well. This means that E°[P(0, id, 0) | -> 0, when t —> oo, 
which implies that 

limP(0,i fc d,0)A fe ^0 (22) 

P°-a.s. for some subsequence {tk : k > 1}, with lim^ tk = oo. Recall that P(0, ifcd, 0) is 
the length of a shortest path from (X(0), 0) (with X(0) = under P°) to {(X(^d), n) : 
n > 0}. Denote one of such shortest paths by <Tk- By the triangle inequality its 
Euclidean length satisfies 

\<Jk\>\0-X(t k d)\>t k -V2s. (23) 

From (22) and (23) one concludes that for any e > and k large enough the length of 
the path o~k is smaller than e time its Euclidean length |<7fc|. Now, (21) follows from 
the fact that the length of the path o~k tends to infinity with fc, which is a consequence 
of tk — » oo and the local finiteness of the graph G (cf Corollary 3.1). 

We conclude the proof by showing that for e small enough, 

^P°{n £ (n)}<oo, (24) 

n 

which by the Borcl-Cantclli lemma implies that n(n) holds P°-a.s. only for a finite 
number of integers n and thus contradicts to (21). To this regard assume constant 
W = w > and let T 5 ™ denote the set of paths a of G of length n, originating from 
(X(0) = 0, 0). Denote also by V§ the analogous set of paths on the graph constructed 
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under assumption W — 0. Note that by monotonicity, 



<pn r -pn 
1 w ' 



(25) 



By the definition 



»°{n £ (n) I $ } = P° (|J{ a G and |a| > n/e } 



(26) 



where the sum bears on all possible n-tuples a = ((Xj 1 , 1), . . . , (Xj n , n)), with Xj i £ 
From this we have 



P°(n £ (n)|$) (27) 

< ^P {<rG:P^,|<r| >n/e|$} 

= E P °{ ff e ^> M > ^ | <7 € VS }P°{ a G 7> " | $} 

< E [Hg ut ' n;W=0 (0)|$]supP°{<7e^,kl >n/e\$,aeVS}, (28) 

where 7ig"*'"' W ' _0 (0) denotes the number of paths of length n originating from (Xq = 
0, 0) under the assumption W = 0. But 

su V P°{a £VZ,\a\>n/e\$,a eVZ} 



< sup E 1 



n 

"[n^j^i-i,™) 



where Jfy = and (i — l,w) = (i — 1) is the indicator of the existence of 

the edge from (Xj i _ 1 , i — 1) to (Xj\ , «) defined by (4), and where we add in the notation 
the dependence on the noise W = w. Using the conditional independence of marks, 
(4), (15) and lack of memory of the exponential distribution of F of parameter (i we 
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have for the path- loss function (1) 



n 



I J >> + '\ jU - 1,0) = 1] 

i=l 
n 

7exp(- M (A|^ i _ 1 -X Ji |)' 3 T«;) . 



i=l 



Hence 



su P P°{a e P.! |cr| > »/e $,cr e } < exp (-fiA^Ti 



where the last inequality follows from a convexity argument. Using this and (27), we 
get 

E°(IL») < E°[H° M *' ,l;W=o (0)]exp(- A/ A /3 Tu;ne-' 9 ) 

< Ccxvi-^Twne- 13 ) 

< exp (»(log(0 - K/e?)) , 

where in the second inequality we used the following result of Corollary 3.1 



Qin,out,n;W=0 



E°[H° 



^•q^j j^out,n;W=0 j^in,k,W=0 ^fe 



and where K is a positive constant. This shows (24) for e small enough, and thus 
concludes the proof. □ 



Proof, [of Corollary 3.5) Without loss of generality assume x — 0. We use the left 
inequality in (6), (14) and the inclusion (16) to obtain 

Px(o),x(y)(0) > L x{0) (0) > T x{0) (0) > f x{0) (0) 

and in consequence 

p(0,y,$)>E[f x(0) (0)|$]. 
Using the isotropy and the strong Markov property of the Poisson p.p. 

E[f A - (o) (0)|<l>]=E o [fo(0)|$y], 
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<&!;§-) is the restriction of $ to the complement of the open ball B = B(q,r) (R), centered 
at (0, R) of radius R > 0, where R is r.v. independent of <& and having for density 

— — 27rAr cxp(— Airr ). 

But since we consider here the SNR graph G 

E o [f o (0)|<%)]>E o [f o (0)|<I>]. 

The result follows now from Lemma 4. □ 



5. SINR space-time graph and routing 

Let us now translate our results regarding the SINR graph into properties of routing 
in ad-hoc networks. 

Firstly, it makes sense to assume that any routing algorithm builds paths on G. 
This takes two key phenomena into account: contention for channel (nodes have to 
wait for some particular time slots to transmit a packet) and collisions (lack of capture 
due to insufficient SINR). 

Our time constant gives bounds on the delays that can be attained in the ad-hoc 
network by any routing algorithms. Of course, realistic routing policies cannot use 
information about future channel conditions. In the case of Poisson p.p. there is hence 
no routing algorithm with a finite time constant. The existence of such an algorithm 
in the case of the Poisson+Grid p.p. remains an open question. In the Poisson p.p. 
case; one can ask about the exact asymptotics of the optimal delay (we know it is not 
linear) and of the delay realizable by some non-anticipating algorithm. 

Let us discuss now the relation of our results to those obtained in [12, 15]. In these 
papers the so called delay-tolerant networks are considered and modeled by a spatial 
SINR or signal-to-noise ratio (SNR) graph with no time dimension. In these models, 
the time constant (defined there as the asymptotic ratio of the graph distance to the 
Euclidean distance) is announced to be finite, even in the pure Poisson case. The 
reason for the different performance of these models lays in the fact that they do not 
take the time required for a successful transmission from a given node in the evaluation 
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of the end-to-end delay. The heavy-tailness of this time (which follows from that of 
the exit time (cf. Proposition 3.1) makes the time constant infinite in the space-time 
Poisson scenario. The reason for the heavy-tailness of the successful transmission time 
is linked to the so called "RESTART" algorithm (see e.g. [17, 18, 1, 16]). In our case 
the spatial irregularities in the ad-hoc network play a role similar to that of the file 
size variability in the RESTART scenario. 
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